Phase separation in the two-dimensional electron liquid in MOSFETs. 



B.Spivak 

Physics Department, University of Washington, Seattle, WA 98195 

We show that the existence of an intermediate phase between the Fermi liquid and the Wigner 
(f) , crystal phases is a generic property of the two-dimensional pure electron liquid in MOSFET's at 

zero temperature. The physical reason for the existence of these phases is a partial separation of the 
uniform phases. We discuss properties of these phases and a possible explanation of experimental 
results on transport properties of low density electron gas in Si MOSFET's. We also argue that in 
certain range of parameters the partial phase separation corresponds to a supersolid phase discussed 
in [17]. 
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<— j ' I. INTRODUCTION 

This work is motivated by experiments [1-13] on transport properties of the two dimensional electron system in high 
mobility Si-MOSFET's at small electron concentration n. These experiments raised doubts about the applicability of 
the Fermi liquid theory and the conventional theory of localization [14,15] to the two dimensional disordered electron 
. ■ liquid at low temperatures. The aim of this article is to prove the existence of zero temperature phases of the two 
dimensional pure electron liquid in MOSFET's which are intermediate between the Fermi liquid and the Wigner 
C| crystal. These phases exist in some interval of concentrations nw < n < n^. The values of the critical concentrations 
I ' nw and til are estimated below. 

This phenomenon is due to a tendency for phase separation which originates from the existence of a first-order 
phase transition between the Fermi liquid and the Wigner crystal phases as a function of n. The difference between 
the crystal-liquid phase transition in MOSFET's and the usual first order phase transitions in neutral systems is 
i i 1 the following. In neutral systems with first order phase transitions the energy of the surface between the phases is 
\ positive and the minimum of the free energy corresponds to a minimal surface area and to a global phase separation. 
In charged systems, like electrons on a positive frozen background, global phase separation does not occur because 
of a large Coulomb energy associated with a non-uniform distribution of electron density. The electron liquid in 
MOSFET's, in a sense, is a system intermediate between these two limiting cases. Similarly to the neutral systems 
with first order phase transitions, the electron liquid in MOSFET's exhibits phase separation. On the other hand the 
iy~j , surface energy of a minority phase droplet of a large enough radius turns out to be negative. As a result at different 
n there is a variety of intermediate phases in this system which are different both from the Fermi liquid and from the 
! Wigner crystal. 

The electron system with phase separation demonstrates a number of unusual features. IfO < (n— nw) *C (til— nw), 
the state of the system corresponds to a small concentration of Fermi liquid droplets embedded into the Wigner crystal. 
The main difference between such a state and the usual Wigner crystal is that it is not pinned by small disorder and 
can bypass obstacles, while the classical crystals at zero temperature are pinned by an infinitcsimally small amount of 
disorder [16]. Phenomenologically this state of matter is similar to the supersolid phase proposed in [17] for the case 
of He 3 and He 4 . The difference is that in our case the origin of droplets of liquid embedded in the crystal is classical 
electrostatic, whereas, in the case [17] the existence of vacancies and interstitials in the ground state of quantum 
O ■ crystals is of quantum origin. 

If < (ni — n) <C {riL — nw), then the state of the system corresponds to a small concentration of Wigner crystal 
droplets embedded into the Fermi liquid. At small concentrations and small temperatures, in principle, these droplets 
can be considered as quasiparticles. 

Droplets of a minority phase interact at large distances via short-range dipole forces rather than via Coulomb forces. 
This means that at T = and at small droplet concentration the system of such "droplet quasiparticles" should be 
in a liquid state similar to He 3 and He 4 which are also liquids at small densities. Thus we can describe the system 
by two-fluid hydrodynamics. However, the statistics of these quasiparticles remains unknown. 

At zero temperature the one-dimensional boundary between the liquid and the solid is a quantum object itself. Due 
to zero-point oscillations of its position there is a region where the wave function has a form which is intermediate 
between the Fermi liquid and the Wigner crystal. Since the electron densities of the Wigner crystal and the Fermi 
liquid are slightly different the fluctuations of the position of the boundary is associated with the fact that the number 
of quasiparticles in the Fermi liquid is not conserved. 
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On the mean-field level this picture of droplet formation in the electron liquid in MOSFET's is similar to the 
partial phase separations which occur in ferromagnetic films [18], charged polymers [19,20], neutron stars [21], doped 
manganites (see for example [22,23]), HTC superconductors [24-27] and 2-dimensional electron systems in the quantum 
Hall regime [28]. All these systems demonstrate a short-ranged tendency for phase separation which is thwarted by 
a long range Coulomb interaction preventing global phase separation. 

The paper is organized as follows. In chapter 2 we show that there is an interval of electron concentrations in which 
the system is unstable with respect to the phase separation. We also estimate the size of minority phase droplets 
embedded into the majority phase and the temperature and magnetic field dependence of the droplet concentration. 
In chapter 3 we discuss transport properties of different nonuniform phases associated with the phase separation. In 
section 4 we compare the theoretical and experimental results on transport properties of the low density electron gas 
in Si MOSFET's. 



II. PHASE SEPARATION NEAR THE POINT OF THE FERMI-LIQUID- WIGNER CRYSTAL PHASE 

TRANSITION 

In this chapter we show that a partial phase separation is a generic property of pure 2-D electron liquids in 
MOSFET's. Consider a two-dimensional electron liquid of density n in a MOSFET separated by a distance d from a 
metallic gate. Electrons interact via Coulomb interaction while a global electric neutrality of the system is enforced by 
the metallic gate with a positive charge density en. The energy density of the system per unit area e(n) = + 
is a sum of the energy density of the capacitor = (en) 2 /2C and the internal energy density of the electron liquid 
e^ el \ In the case of a uniform electron distribution the capacity per unit area is C = Co = 1/d. 

At high electron densities na 2 B 3> 1 the kinetic energy of electrons is larger than the potential energy and the 
interaction can be taken into account by a perturbation theory. (Here as is the electron Bohr radius.) In this case the 
system can be described by Fermi liquid theory, the difference between the effective m* and the bare m electron masses 
is small, and = ejf 1 ~ n 2 /m. On the other hand, in the opposite limit na 2 B <C 1 (but still nd 2 ^> 1) the potential 
Coulomb energy of electrons is much larger than the kinetic energy and the ground state of the system is a Wigner 

crystal with = = —e 2 r?l 2 (see, for example, [29]). Thus at zero temperature there is a critical electron 
concentration n c where the phase transition between the Fermi liquid and the Wigner crystal phases takes place. 
According to Landau mean field theory this transition is of the first order (see for example [30]). The n-dependence 
of the energy densities of the Fermi liquid ez,(n) = e c + and the Wigner crystal ew{n) — + effl P nases near 
the critical density n c is shown schematically in Fig. la. 

In the limit of small densities nd 2 <C 1, due to the existence of the image charges in the gate, the interaction 
between adjacent electrons has a dipole character. In this case the ratio between the potential and the kinetic energy 
decreases as n decreases. Therefore, the small electron n the electron system is a weakly interacting Fermi liquid. 
Thus we arrive at the conclusion that there exists another critical point r^P <~ 1/d 2 which corresponds to a second 
Wigner crystal- Fermi liquid transition. The phase diagram of the electron system at T = is shown in Fig. lb. If 
d < d* ~ 38as, than the system is in the liquid state at any value of n. Here the factor 38 is the result of numerical 
simulations [31]. 



A. The mean field description of the phase separation. 

In the approximation when C = Co the qualitative picture of the phase transition is the same as the picture of any 
first order phase transition in neutral systems. In particular, there is an interval of electron densities % < n < Ul 
shown in the Fig. la where there is a phase separation, which means that there is a spatially nonuniform distribution of 
the Wigner crystal and Fermi liquid phases coexisting in equilibrium. In the case of large d one can linearize ef^w (n) 
near the point n = n c . As a result, we have 

, (p-w - Ml) x 
n L ,w = n c ± — — (1) 

where fi W x = {de^ L /dn)\ n= „ c . 

One can get from Eq.l an estimate n c as/ d for the size of the interval of electron densities where the phase separation 
occurs. Values of d/as in various MOSFET's range from of order one to 50. 
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The relative fractions of these phases Xw and xl are determined by the Maxwell rule. At (jil — n) <C (til — nw) 
the fraction of the area occupied by the Wigner crystal xw "C 1 is small while in the case (n — nw) <C (til — nw) the 
fraction of the area occupied by the Fermi liquid xl <C 1 is small. 

n-n w .L 

x w ,L = ± (2) 

n c 

The compressibility of the system v = d 2 e/dn 2 should exhibit jumps of order e 2 d at points n — til, nw- 
The crucial difference between first order phase transitions in neutral systems and in the system of electrons in 
MOSFET's arises when one considers shapes of the minority phases. In the case of neutral systems the surface energy 
density a is positive. Therefore in equilibrium the system should have a minimal area of the surface separating the 
phases, leading to global phase separation. On the other hand, in the three dimensional charged systems the global 
phase separation is impossible because of the large Coulomb energy associated with the charge separation. It is 
possible, however, that in this case the electron system consists of bubbles and stripes of different electron density 
[24,27], provided the tendency for phase separation is strong enough. 

The situation in MOSFET's is very different. On one hand, in the approximation when C — Co global phase 
separation is possible at an arbitrary value of (/iw — Hl)- On the other hand, it turns out that for large droplets of 
the minority phase the surface energy is negative. To prove this one has to take into account the finite size corrections 
to the standard formula for the capacitance [32] 

_ „ R, WttR 

C = C + -In — (3) 

where A and d are the capacitor area and thickness respectively and R = \f~A is the capacitor size. Consider, for 
example, the case when xw <C 1- Then xw can be determined by the Maxwell rule in the approximation when the 
second term in Eq.3 is neglected and C = Co- Expanding e(n) with respect to the second term in Eq.3 and taking 
into account also the microscopic surface energy we have an expression for the energy of the surface 



E( SU rf) = -7^N w e 2 {n w - n L ) 2 d 2 R w In 1&7r ^ w + N w <t2ttR w (4) 



We assume that the Wigner crystal phase embedded into the liquid consists of droplets of radius Rw and concentration 
Nw and take into account that inside the droplet n ~ nw- Thus, at large Rw the surface energy Eq.4 turns out to 
be negative. We have to find a minimum of Eq.4 at a given total area occupied by the minority phase, which gives us 
the characteristic size of the droplet 

Rw ~ 4c el ( 5 ) 

107T 

with 7 = (e 2 a)/2n(fiw — ^l) 2 - A similar expression was obtained in [33] for a different problem. 

The analogous calculation for the case xl -C 1 gives the expression for the radius of liquid droplets embedded into 
the crystal which is identical to Eq.5. 

At the point of the transition the values of a and {a w — ul) 2 /e 2 arc of the same order and at present nothing is 
known about the value of the dimensionless parameter 7. Even the fact that a > is not proven. I would like to 
also note that in the case of the first order phase transitions which are close to the second order one we always have 
7 < 1. 

In this article we assume that 7 > 1. To illustrate the physical meaning of this inequality we consider the case 
when 2-D electron liquid is compensated by a uniformly charged positive frozen background with a charge density 
en. In this case the Coulomb energy of a droplet associated with the phase separation is, roughly, R/d times larger 
than in the MOSFET's case. The most dangerous point with respect to the phase separation instability is n = n c 
(see Fig. la). For example, let us compare the energies of the uniform liquid state with n = n c and a nonuniform 
state which contains two droplets embedded into the liquid. The first droplet is a liquid with electron concentration 
n\ = n c + Sn, while the second term is a crystal with electron concentration n2 — n c — Sn. Suppose the droplets have 
the same radius R. Linearizing e.L,w{n) with respect to Sn we estimate the energy difference SE between these two 
states as 

SE ~ {p L - n w )TrR 2 Sn + ( eSn7TR2 } 2 + 2txRg (6) 

R 

The first term in Eq.6 corresponds to a decrease of the energy due to the phase separation. The second one corresponds 
to the positive Coulomb energy associated with the nonuniform distribution of the electron density and the third term 
is the surface energy. A minimization of Eq.6 with respect to Sn gives us Sn <~ (u w — n L )jRe 2 and 
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SE mm ~(2na- [fXW / L> )R (7) 

The assumption 7 > 1 means that E m i n in Eq.7 is positive and that 2-d electron liquid on a, frozen positive background 
does not exhibit a phase separation. (We would like to metion that if the value of \nw — Ml| is big enough a microscopic 
phase separation in the charged liquid on frozen positive background can take place even in 3D case. Such situation 
has been considered in [24-27,58]) 

On the mean field level our problem is similar to [18,19]. Using this analogy we conclude that in the middle of the 
interval (nw,nL) there is a stripe phase. The phase diagram of the system is shown schematically in Fig. 2. The main 
difference with [18,19] is the following. In [18,19] all phase transitions between uniform, bubble and stripe phases are 
of the first order, whereas in our case the transitions between uniform (Fermi liquid and Wigner crystal) phases and 
the bubble phases are continuous. The transitions between the bubble phases and the stripe phase would be the first 
order one. However, such a transition would have an interval of concentrations where phase separation would take 
place. In this case the presented above arguments could be repeated. Thus we expect more complicated structures 
than bubbles and stripes phase to exist between the bubble and the stripe phases. Since the complete solution of this 
problem remains to be found we indicated this in Fig. 2 by shaded lines. 

Let us now estimate dependence of Xw,l{T, Hu) on the temperature T and the magnetic field Hu parallel to the 
film. It is determined by the corresponding dependence of the free energies for the Fermi liquid and Wigner crystal 
phases. At small T and Hu one can neglect the T and Hu dependences of ew,L and we have the following expression 
for the free energy densities of the liquid and the Wigner crystal phases 

FwA H \\) = e wx - M w , L H l{ n - TnS WtL (8) 

where Sw and Sl are the entropies of the crystal and the liquid phases respectively, while Mw and Ml are the 
corresponding spin magnetizations per electron. 

As a result, one can obtain how xw,l(T, H\\), and n c w,L (T, Hn ) depend on T and H\\ by making the following 
substitution in Eqs.1,2 

(Hw - Ml) -» (Pw - Ml) - (M w - M L )H\\ - T(S W - S L ) (9) 

At small HbHu «T« Ep we have Mw,l = Xw,lHu, where xw and \l are linear susceptibilities of the crystal and 
the liquid respectively. (Here /ib is the Bohr magneton and Ep is the Fermi energy.) At low temperature T <C Ep 
the spin susceptibility of the Wigner crystal xw ~ Mb/^" ^ XL is much larger than the spin susceptibility of the 
Fermi liquid. The entropy of the crystal Sw ~ In 2 » Sl ~ T/Ep is mainly due to the spin degrees of freedom and 
much larger than the entropy of the Fermi liquid. Thus xw increases linearly with T and quadratically with Hu , 
which means that both the temperature and the magnetic field parallel to the film drive the electron system toward 
the crystallization [34]. (We assumed that the temperature is larger than the exchange energy between spins in the 
Wigner crystal). These effects are known in the physics of He 3 as Pomeranchuk effects. 

In the intermediate interval of magnetic fields T < ^bH\\ < Ep spins in the Wigner crystal are completely polarized 
while the Fermi liquid is still in the linear regime. In this case xw increases linearly with Hn . 

At high magnetic field H\\ > Hjr <~ Ep/fiB both Fermi liquid and Wigner crystal are spin polarized and xw(T, H\\) 
saturates as a function of Hu. We assume that £l(H\\ = 0) < el(H > H^) and, therefore Xw{H\\ = 0) < xw{H > 
On the other hand, the spin entropy of the Wigner crystal is frozen in this case. As a result, at Hu > H^ the 
temperature dependence of Xw,l(T, Hu) is suppressed significantly. 

B. Quantum properties of the droplets of minority phase embedded into the majority one. 

In principle, at small enough concentrations and at small temperatures droplets of the minority phase embedded 
into the majority one should behave as quasiparticles. Since the system is translationally invariant, they should be 
characterized by momentum (or by quasi- momentum) . The momentum coincides with the flux of mass. Thus these 
quasiparticles carry a mass M*, a charge eM*/m and a spin. The characteristic temperature of quantum degeneracy 
\sT* ~ N w /M*. 

The value of M* depends on the mechanism of motion of the droplets, which in turn depends on whether the 
surface between the crystal and the liquid is rough or smooth. 

Consider for example the case of Wigner crystal droplet embedded into the liquid. In the case of a smooth surface, 
motion of the droplet is associated with a redistribution of the liquid mass on the distance of order Rw- In this case 
we can estimate the effective mass of the droplet as 



4 



M* ~ mn c nR^ 



■w 



(10) 



In the case of rough surfaces the motion of the droplet is associated with melting and crystallization of different parts 
of it. Since the (ul — nw) <C nc the liquid mass to be distributed and, consequently, effective mass of the droplet 



in this case is much smaller than Eq.10. 

Droplets of the minority phase interact at large distances via short-range dipole forces rather than via Coulomb 
forces. At small enough concentration of the droplets the amplitude of quantum (or classical) fluctuations of their 
positions is larger than the typical distance between them. Thus the liquid droplets are distributed uniformly over 
the whole crystal. In other words, at T = the system of such "droplet quasiparticles" should be in a liquid state 
similar to He 3 and He 4 which are also liquids at small densities. Thus we can describe the system by two-fluid 
hydrodynamics. In this case statistics of the "droplet quasiparticles" becomes important. In this respect we would 
like to mention a difference between the droplets of the liquid embedded into the crystal and the droplets of the 
Wigner crystal embedded into the liquid. 

a. The droplets of the liquid are topological objects which, in principle, arc not different from vacancies or inter- 
stitials in quantum crystals He 3 and He*. In order to create such objects in Wigner crystal one has to add or to 
remove from the lattice an integer number of electrons. Therefore, the liquid droplets have a definite statistics: they 
are either fermions or bosons [17]. 

The main feature of the phase where there are droplets of liquid embedded into the crystal (supersolid) is it's 
ability to bypass static obstacles. In other words, unlike conventional crystals supersolids are not pinned by disordered 
potential of small amplitude. This will manifest itself in the finite conductivity of the system. 

From the phenomenological point of view this is very similar to the scenario of "supersolid" which has been 
introduced by A.F.Andreev and I.M.Lifshitz [17] for quantum crystals of helium near the quantum melting point. 
They assumed that the crystals contain zero point defects (vacancies or interstitials) in the ground state and therefore 
the number of atoms and number of sites in the crystals are different. The difference with [17] is that the origin of 
the negative surface energy Eq.4 is purely classical. Conversely, following [17] the existence of point defects in the 
ground state could be of quantum origin. Namely, the kinetic energy of the point defects can be larger than the 
energy required for their creation. Thus, the supersolid phase [17] can be considered as a particular case of a more 
general situation of the phase separation when the radius of liquid droplets embedded into the crystal is of order n~ 2 . 
This would mean that the surface energy is renormalized to a small (or negative) value. Indications of the existence 
of such a phase have been reported in numerical simulations [35] . 

b. The case of droplets of the Wigner crystal embedded into the Fermi liquid is different because they arc not 
topological objects. In principle such droplets can contain an additional charge and spin which can be fractional or 
even irrational. A fundamental problem associated with this fact is that statistics of such quasiparticles is unknown. 

To illustrate this point we consider a process of tunneling between two states: a state of uniform Fermi liquid and 
a state when there is one crystalline droplet embedded into the Fermi liquid. These two states have different total 
electron charge. Thus the tunneling between these states is associated with a redistribution of this charge to (and 
from) the infinity. It is important that the action S associated with this process in the pure two dimensional case is 
finite. One can estimate it in a way similar to [36]. On distances larger than the droplet size R\y one can write the 
action in terms of the time-dependent electron density n(r, t) 



Here t is the imaginary time and M* is given by Eq.ll. We approximate that n(r, t) <~ M* /mR 2 (t) at |r| < R(t) and 
n = at |r| > R(t). Eq.12 corresponds to the potential energy contribution to the action. As usual, the contribution 
from the kinetic energy is of the same order. Assuming that R(t) — vpt we get an estimate 5* ~ (eM*) 2 d/m 2 RwVF- 
Thus, in principle, the wave function of the object is a coherent superposition of the wave functions of a uniform 
Fermi liquid and a Wigner crystal droplet. In this situation it is quite likely that the additional charge associated 
with such an object is not an integer. This is the reason why the nature of the ground state of the system remains 
unknown. 

The quantum melting of the phases, which are intermediate between the bubble and the stripe phases is even more 
complicated and we leave this question for further investigation. 



M* ~ m(nL — nw)^R%, 



■w 



(11) 




(12) 
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III. TRANSPORT PROPERTIES OF THE ELECTRON SYSTEM WITH DROPLETS OF A MINORITY 

PHASE EMBEDDED INTO THE MAJORITY ONE. 

In this section I will consider cases when quantum statistics of the system of droplets of the minority phase is not 
important. 

The electron-electron scattering conserves the total momentum of the electron system and therefore does not 
contribute to the resistance of the system. To estimate it we have to consider the electron system in the presence of 
a random elastically scattering potential. 

The electron transport picture in the electron liquid with partial phase separation is quite rich. In particular, there 
is a region of electron concentrations where the hydrodynamics of the electron liquid is similar to the hydrodynamics 
of the liquid crystals [26]. In this paper we consider only cases where either there are crystalline droplets of small 
concentration embedded in the liquid (xl <C 1), or there are liquid droplets with xw *C 1 embedded into the crystal. 
In these situations, in principle, there are two types of current carriers in the system: electron quasiparticles and 
charged droplets of the minority phase. In this article we will ignore the contribution of the droplet motion to the 
charge transport. To illustrate the possible T and H dependence of the resistance we consider below only several 
limiting cases leaving a detailed analysis for future investigation. 



A. The case when crystal droplets of small concentration are embedded in the electron liquid. 

Let us consider the case xw <C 1 when crystalline droplets of small concentration are embedded into a Fermi liquid. 
We will assume here that the Wigner crystal droplets are either pinned by a small scattering potential, or have a 
short mean free path. We also assume that otherwise the impurities do not affect the thermodynamic properties of 
the system. The contribution to the resistance of the system from the scattering of quasiparticles on droplets has the 
form 

where kp is the Fermi momentum of the Fermi liquid, l( e ,w) = 1/NwRw is the quasiparticle mean free path, and 
Nw — xw I Rw i s the concentration of droplets of the Wigner crystal. Thus as followes from Eq. 1,2, 9, 13 at small T 
the resistance of the electron system increases linearly in T. At small H\\ it increases quadratically in H\\, while in 
the intermediate interval of H\\ it increases linearly in H\\ . The saturation of the magnetoresistance as a function of 
takes place at ify > HZ when the electron Fermi liquid gets polarized. 

At H\\ > ify the spin entropy of the Wigner crystal is frozen. Therefore, as it has been discussed, xw{T) and the 
resistance of the system do not have a significant T-dependence. 

The H\\ dependence of the resistance p(Hn) of the metallic phase at small T is shown schematically in Fig. 3a. The 
T dependences of p(T) at H\\ = and H\\ > £Z| are shown in Fig. 3b. 

Eventually at high enough temperatures the crystalline droplets melt. Since at this point r s » 1 the melting 
temperature T m <C f2 p is much smaller than the plasma frequency at the wave vector of order of the inverse inter 
electron distance. Here rs is the ratio between the potential and the kinetic energies of electrons. Let us now discuss 
the T-dependence of p(T) in this temperature interval. Though in this case the liquid is not degenerate, it is strongly 
correlated. Therefore the electron-electron scattering in the liquid is very effective and the local equilibrium is reached 
in a short time on a spatial scale of order n -1 / 2 . As a result, the flow of the electron liquid near an impurity can 
be considered in the framework of hydrodynamics. In the two-dimensional case the moving electron liquid exerts a 
force on an impurity, which is given by the Stokes formula F ~ r\uj In ((77 / 'nua)) , [37]. Here u, r) and a are the liquid 
hydrodynamic velocity, viscosity of the electron liquid and the impurity radius respectively. In a system with a finite 

1/2 

concentration of impurities the logarithmic factor in the equation for F should be substituted for In(l/oJV f ), where 
Ni is the concentration of impurities. Thus the resistance of the electron system has the form [39,38] 

The viscosity of the strongly correlated liquid in the semi-quantum regime has been considered theoretically in [40] 
for the case of liquid He 3 . It was conjectured that 

(15) 
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We can apply this result to the case of electron liquid in the semi quantum regime (T m <T< Q p ) as well. Thus we 
arrive at the conclusion that at high temperatures the resistance should decrease inversely proportional to T and that 
it should have a maximum at T ~ T m ~ Ep. It is interesting to note that, as far as I know, the experimental data 
on the T-dependence of the viscosity of He 3 in this relatively high temperature region are unavailable. However, we 
can look at data for the viscosity of He 4 , which in this temperature interval is supposed to be similar to He 3 [40]. 
Though the experimental data for He 4 are in a reasonable agreement with Eq.15, we would like to mention that the 
viscosity of He 4 changes only by a factor of two in the temperature interval between Ep and the evaporation point. 

B. Strongly correlated Fermi liquid in the presence of a scattering potential. 

Let us consider the case n > nw, but r s ^> 1. Then at Ep -C Q p the electron system is a strongly correlated 
Fermi liquid. The main feature of such a liquid is that at small distances and at small (imaginary) times it behaves 
like a solid [34]. It has been suggested in [41] that the cross-section of quasiparticle scattering on a short ranged 
impurity with a radius of order n~ 2 , is significantly enhanced by the electron-electron interaction. The nature of 
the enhancement becomes especially clear if we consider the interval of electron densities close to the critical point 
< n — riL <C rip and the case when the fluctuations of the external potential have a relatively small amplitude. 
Then the system can become split into the regions of a Fermi liquid and a Wigner crystal. I would like to mention 
that the linear T and H\\ dependences of the resistance mentioned above are generic for strongly correlated electron 
system and are valid in this case as well. 

The fractions of volume occupied by the Fermi liquid and the Wigner crystal depend on n and therefore the system 
should exhibit a percolation-type zero-temperature metal-insulator transition as n decreases and the area occupied by 
the Wigner crystal grows. There is however a significant difference with respect to the percolation transition, which 
originates from the fact that the position of the Wigner crystal-Fermi liquid boundary has quantum fluctuations. It 
is these fluctuations which determine the character of the electron transport near the transition point. The properties 
of the current carriers in this region are very different from properties of the Fermi liquid quasiparticles. 

The significance of quantum fluctuations becomes even more clear if we consider quantum properties of the surface 
between the Wigner crystal and the Fermi liquid. At zero temperature the surface is a quantum object itself. There 
are at least two scenarios for the state of the surface at T — 0: it could be quantum smooth or quantum rough. In 
the first case the excitations of the surface are essentially the Ralaigh surface waves which conserve the charge inside 
the droplets. In the second case there is a new type of excitations at the surface: crystallization waves which do not 
conserve the total charge inside the droplets [42]. 

The problem of quantum roughening has been discussed in the framework of the properties of the boundary 
between solid and liquid He 3 and He 4 [42]. In the case of 2-d surfaces between a 3-d quantum liquid and crystal it 
has been argued that the surface is always quantum smooth [43]. At the moment nothing is known about the state 
of the boundary between 2-dimensional liquid and solid. We would like to mention, however, that quantum effects 
(including the quantum roughening) arc more pronounced in the case considered above because it is two-dimensional 
and because the jump of the electron density in this case is small. In any case, due to quantum fluctuations, there is 
a region near the boundary whose properties are intermediate between the liquid and the solid properties. 

In the conclusion of this section we would like to mention that the linear in T increase of the resistance at small T 
is a generic property of the model. 

IV. A COMPARISON BETWEEN PREDICTIONS OF THE THEORY AND EXPERIMENTAL RESULTS 

IN SI MOSFET'S. 

A. A review of experimental results on Si-MOSFET's. 

In this subsection we present a short list of experimental results on the high-mobility two-dimensional electron 
liquid in Si MOSFET's [1—13] which seem to be in contradiction with Fermi liquid theory and with the conventional 
single particle localization theory of disordered 2-dimensional conductors [14,15]. 

A. The electron system exhibits a "transition" as a function of n from a metallic phase, where the resistance of the 
system saturates at low temperatures, to an insulating phase, where the resistance increases as T decreases. The value 
of the critical concentration ni MI ^ of the transition depends on the amount of disorder in the sample and corresponds 
to r s = » 1. Here r s is the ratio between the electron potential and kinetic energies. 
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B. At T = and for the electron concentration sufficiently close to the critical one, increasing the magnetic field 
parallel to the film drives the system toward the insulating phase [3,5,12]. Thus the critical metal-insulator 

concentration n c (H») increases with H\\. 

In the metallic phase (n > n M1 (Hn =0)) and at small T the system exhibits a big positive magnetoresistance as a 
function of H||. This magnetoresistance saturates at ffy > H^(n) and p(H^)/p(0) 3> 1 [3,5]. 

C. In the metallic phase at Hn = and T < Ep the resistance p(T) significantly increases with increasing temper- 
ature. The characteristic value of <iln p/dT > Ep 1 at small T is large and depends on the value of n — n c . 

D. If at if || > if| the system is still in the metallic phase (n > n c (if||)), the T-dependence of the resistance is much 
weaker than in the ff|| = case [6,8,9]. 

E. The value of fi| decreases significantly as n approaches 

B. Qualitative explanation of experimental results. 

In this subsection we present a qualitative explanation of the experimental results [1-13] 

A. The existence of the metal-insulator phase transition. 

The theoretical picture presented above involves a transition between the liquid and the crystal as a function of 
n. Therefore it can explain qualitatively the existence of the metal-insulator transition observed in the experiments. 
Namely, the fractions of volume occupied by the Fermi liquid and the Wigner crystal depend on n and therefore 
the system should exhibit a percolation-type zero-temperature metal-insulator transition as n decreases and the area 
occupied by the Wigner crystal grows. The transition takes place when the Wigner crystal droplets overlap and block 
the electron transport through the Fermi liquid area. 

The experimental values of ni^ 11 - 1 correspond to r s <~ 10 — 20. At present it is difficult to say how close this 
value is to til, or n c in the pure case. The critical value for the transition r SiC = 38 [31] was obtained by numerical 
simulations. However it can not be applied to the case of electrons in Si MOSFET's because of the existence of two 
almost degenerate electron valleys. Another reason for possible inapplicability of the results of [31] to Si MOSFET's 
is that the calculations [31] were restricted to the case of zero temperature while the experiments have been performed 
at temperatures larger than the spin exchange energy in the Wigner crystal. Thus the Pomeranchuck effect has not 
been taken into account in [31]. Finally in the critical value of r s can be different in the disordered case. 

B. The positive magnetoresistance of the metallic phase in the magnetic field parallel to the film. 

The large positive magnetoresistance of the metallic phase in the parallel magnetic field is connected to the fact that 
Xw ^> Xl and therefore the magnetic field parallel to the film drives the electron system toward the crystallization 
[34]. (See Eqs.1,2,9,13). The magnetoresistance should saturate when ff|| > if| and the electron Fermi liquid is 
polarized. 

C. The temperature dependence of the resistance in the metallic phase. 

The significant increase of the resistance as a function of temperature can be explained naturally as a consequence 
of the Pomeranchuk effect: The spin entropy of the Wigner crystal is larger than the entropy of the Fermi liquid and, 
therefore, the Wigner crystal regions grow with increasing temperature. 

At high temperatures the droplets of crystal melt. It follows from Eqs.13,15 that in this temperature range the 
resistance decreases with increasing T. It is unclear at present whether the experiments support this picture. 

D. The temperature dependence of the resistance in the metallic phase at large Hn. 

The Pomeranchuk effect disappears when ff|| > HS and electron spins are fully polarized. In this case entropies of 
both the liquid and the solid are much smaller than the spin entropy of the crystal at ff|| = 0. This means that in 
the leading approximation the areas occupied by the crystal and the liquid are T-independent. This explains the fact 
that in the metallic state at H\\ > iff? the T-dependence of the resistance is much smaller than in the case ff|| = 

[8,9]. (The ratio dp/dT(H\\ = Q)/dp/dT(H\\ > £f|) can be as big as 10 2 ). 

E. The n-dependence of H?,. 

Perhaps the most direct check of the concept of the Fermi liquid which is close to crystallization is the measurement 
of the n-dependence of the magnetic field H^(n) which polarizes the liquid. In the case of a non-interacting Fermi 

liquid ff| = i2| (0) = E f /hb is a smooth function of n. The problem of the n dependence of the critical magnetic 
field H?,(n) in strongly correlated liquids near the crystallization point and the origin of the strong enhancement of 

the spin susceptibility has been discussed in the context of the theory of liquid He 3 [34] . It has been pointed out that 
there are two different scenarios for the origin of the significant (factor 15) enhancement of the spin susceptibility of 
He 3 near the crystallization point. 
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a) The system is nearly ferromagnetic which means that it is close to the Stoner instability. In this case the linear 
spin susceptibility \l is large, but other coefficients a m in the expansion of the energy 

e L = X L 1 M 2 + a i M i + ..a m M m (16) 

with respect to the spin magnetization M are not small. Here m is an even integer. In this case fl| ~ , which 
is relatively large. 

b) The system is nearly solid. In this case both Xi} an d other coefficients a m in the Eq.15 decrease significantly as 
n approaches the crystallization point n c . In this case i2|(n) < i2| (0) is small. 

In the case of He 3 the value of H?An) has never been measured. In the case of electrons in Si-MOSFET's it has 

been measured in [13]. A dramatical decrease of H^(n) compared to H^ ^ has been observed as n approaches n c . In 
our opinion these experimental results support the model of a nearly solid Fermi liquid which is at r s 3> 1. Conversely 
it is unlikely that the system is close to the Stoner instability. 

C. A comparison with alternative explanations of experiments on transport properties of the metallic phase 

of the electron system in Si MOSFET's. 

In this section we compare the explanation presented above with another explanation given in [44-49] . It is based 
on the fact that a single short range impurity in a metal creates Friedel oscillations of the electron density. Due to the 
electron-electron interaction the quasiparticles in the metal are scattered not only from the impurity but also from 
the modulations of the electron density. At finite temperature the Friedel oscillations decay exponentially at distances 
larger than the coherence length of the normal metal vf/T. As a result, at low temperatures (p(T) — p(0)) <~ CT with 
C > [44-47]. The exchange contribution to the resistance has not been taken into account in [44-47]. It has been 
shown [49] that in the presence of the exchange interaction at t _1 <C T <C Ep the quantity (p(T) — p(0)) remains 
linear in T. However, at r s <C 1 the coefficient C < has negative sign, which is different from [44-47]. On the other 
hand, the expriments were performed in the regime r s > 1. They yeald a positive value of the coefficient C > 0. At 
finite value of r s ~ 1 the theory [49] predicts, that the coefficient C chenges it's sign agein and become positive. 

At relatively high temperatures T ~ Ep one can neglect the interference corrections and the temperature dependence 
p(T) is determined by the corresponding dependence of the thermal velocity and the electron scattering cross-section 
in a nondegenerate gas. In this regime p(T) decreases with increasing T [44,47]. At this point we would like to 
mention that at low T Eqs. 1,2,9, 13 also predict the increase of the resistance linear in T as well. Eqs.14,15 also 
predict the existence of the maximum of p(T) at T ~ Ep and decrease of the resistance at T > Ep. Thus, both the 
theory presented above and [44-49], in principle, could explain qualitatively the T-dependence of the resistance of 
the metallic state. 

The situation with the magnetoresistance in the parallel magnetic field is more delicate. Strictly speaking, the 
interference corrections to the Drude conductivity calculated in [44 49] are relevant only at r s <C 1 and at small 
T and H\\, when the effects are small. On the other hand, at H\\ > H$, when the effects are large, the interference 
corrections are irrelevant, and the value of the magnetoresistance (p(H») — p(0)) is determined by the Hu dependence 
of the Drude part of the resistance, which is due to the H\\ dependences of the Fermi momentum and the scattering 
cross-section of quasiparticles. (This part of the magneto resistance has not been taken into account in [49]). In this 
case a single electron theory yields a big and negative Drude magnetoresistance in contradiction with the experimental 
fact that it is big and positive. 

In connection with this I would like to make several points. 

a. At r s > 1, the diagrammatic calculations are not under control. In the framework of the conventional diagram 
technique it is difficult to account for all effects associated with the strong correlations at r s 3> 1, including the giant 
renormalization of the electron scattering cross-section on impurities, the effects of phase separation, existence of the 
crystallization waves at the boundary between the two phases, and, finally, the Wigner crystallization itself. 

b. The mechanism considered in [44-49] can not explain the increase of the resistance as a function of Hu and 
T which is significantly larger than unity. This is because the amplitude of the potential created by the Friedel 
oscillations of the density created by an impurity potential is smaller than the impurity potential itself. This theory 
also can not explain why the temperature dependence of the resistance p(T) is suppressed so dramatically by the 
magnetic field parallel to the film. On the other hand, the theory presented in this article can explain these facts. 

c. These mechanisms of T and H\\ dependence of the resistance are based on very different physics. This can be 
seen, for example, from the fact that all single electron interference phenomena including the Friedel oscillations are 
smeared by finite temperature. Conversely, the fraction of the Wigner crystal increases with temperature. 
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d. The amplitude of the Friedel oscillations is suppressed significantly in the case when the scattering potential is 
a smooth function of coordinates on the scale of the electron wavelength. This is exactly what happens when the 
scattering cross-section is significantly renormalized by the fact that near a short range impurity there are crystalline 
droplets and the position of the crystalline surface exhibits quantum fluctuations. Thus, in a sense, the mechanisms 
based on single electron interference and the mechanism based on quantum fluctuations of the solid-liquid boundary 
compete with each other. 

In order to distinguish between these two mechanisms one needs to perform experiments on samples with higher 
mobility, where effects considered in this article will be much larger than unity. 

Finally we would like to mention that the theory [44-49] may be relevant to experiments on the two dimensional 
electron system in GaAS samples [50,51]. 



V. CONCLUSION 



We have shown that due to the existence of metallic gates in MOSFET's the phase separation is a generic property of 
pure electron liquids. The proof is based only on the assumption about the existence of the first order phase transition 
between the uniform Fermi liquid and the Wigner crystal phases and on electrostatic properties of two dimensional 
electron system. This distinguishes the theory presented above from the theories [52,38], which attempted to explain 
the experiments using the fact that in the 2-d electron liquid there is a first order phase transition between the 
Fermi-liquid and the Wigner crystal which is destroyed by small disorder. This difference, however, manifests itself 
only at relatively small values of d and at relatively small amplitude of the disorder. Qualitative pictures of the T 
and independences of the resistance of the " metallic" phase are, roughly speaking, the same for the model presented 
above and for those considered in [38]. 

It is an open question how a disorder of finite amplitude affects the results presented above. In some regimes 
the system can demonstrate a glassy behavior characteristic for crystals in the presence of disorder. Experimental 
indications of glassy behavior of the electronic system in Si MOSFET's have been reported in [53,54]. 

In this paper we considered only bubble phases which exist near the critical concentrations and nw- In the 
interval til < n < nw the system, will probably exhibit a sequence of quantum phase transitions. In particular, it is 
likely that at electron densities close to n c there is a stripe phase, which is similar to [24,27,28]. 

In conclusion we would like to mention that the picture presented above is in many respects similar to the quantum 
critical point of strongly correlated electron systems considered in [55-57]. In particular, the Fermi liquid state with 
densities close to nw will demonstrate very large sensitivity to imperfections, which is characteristic for the " almost 
critical" quantum state [55]. 

This work was supported by Division of Material Sciences, U.S. National Science Foundation under Contract No. 
DMR-9970999. We would like to thank A.F. Andreev, E. Abrahams, S.Chakrovarty, A.Efros, M.Gershenson, S. Kivel- 
son, S. Kravchenko, L. Levitov, D. Maslov, V. Pudalov, B. Shklovskii, M. Sarachik, and S. Vitkalov for useful discus- 
sions. 
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a) b) 
FIG. 1. a. The dependence of the energy densities of the Wigner crystal and the Fermi liquid phases ew,L(n) on the electron 
density n. Symbols W and L correspond to the Wigner crystal and the Fermi liquid phases respectively, b. The effective phase 
diagram of the 2D electron system at zero temperature. 
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FIG. 2. The phase diagram of the 2D electron system at T = 0. Symbols WC and FL correspond to the Wigner crystal and 
the Fermi liquid phases respectively. The shaded regions correspond to phases which are more complicated than the bubble 
and the stripe phases. 
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FIG. 3. a)The H\\ dependence of the resistance p{H\\). b) The temperature dependence of the resistance p(T). 
line corresponds to the case H» = 0, while the dashed line corresponds to the case H\\ > H?, 
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